Multipole Expansion for Relativistic Coulomb Excitation 



H. Esbensen 

Physics Division, Argonne National Laboratory, Argonne, Illinois 60439 

C. A. Bertulani 

Physics Department, Brookhaven National Laboratory, Upton, NY 11973-5000 - USA 

(February 8, 2008) 

Abstract 

We derive a general expression for the multipole expansion of the electro- 
magnetic interaction in relativistic heavy-ion collisions, which can be em- 
ployed in higher-order dynamical calculations of Coulomb excitation. The 
interaction has diagonal as well as off-diagonal multipole components, asso- 
ciated with the intrinsic and relative coordinates of projectile and target. A 
simple truncation in the off-diagonal components gives excellent results in 
first-order perturbation theory for distant collisions and for beam energies up 

to 200 MeV/nucleon. 

PACS number(s): 25.70.-z, 25.70.De 
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I. INTRODUCTION 



The Coulomb excitation of nuclei at relativistic energies is commonly calculated in first- 
order perturbation theory, using the formalism developed in Ref. The derivation is 
restricted to distant collision, where the intrinsic radial coordinates are smaller than the 
minimum distance between the colliding nuclei. In calculations of the Coulomb dissociation 
of proton halo nuclei, it is of interest also to consider close collisions, since the density of the 
valence proton can extend to very large distances. Moreover, higher-order processes may 
also play a role, as suggested by the non-relativistic calculations of the breakup reported 
in Ref. §. 

In order to calculate the Coulomb dissociation of halo nuclei at energies where relativistic 
and higher-order effects cannot be ignored, one would need realistic multipole form factors for 
the electro-magnetic interaction. Such form factors have recently been studied but they 
were only determined in the long-wave-length limit for distant collisions. We therefore find it 
timely to study the multipole decomposition of the interaction at relativistic energies, both 
for close and distant collisions. The form factors we derive can be applied in higher-order 
dynamical calculations of two-body breakup reactions, such as the continuum discretized 
coupled-channels calculations Q], or in numerical methods that evolve the two-body wave 
function in coordinate space 

We base our study on the so-called Lienard-Wiechert potential and include the effect 
of the convection current. The magnetization current, on the other hand, is ignored for 
simplicity. We assume that the relative velocity of the interacting nuclei is a constant, so 
that the retardation effect associated with a time- dependent velocity can be ignored. 

The multipole expansion of the interaction contains both diagonal and off-diagonal mul- 
tipole components. The off-diagonal components appear because of the Lorentz contraction. 
Our approach is similar to the method developed by Baltz et al. in their study of the 
e~^e~ pair production in relativistic heavy-ion collisions. We extend their study and derive 
the multipole fields that are relevant to Coulomb excitation. 
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The multipole expansion of the Lienard-Wiechert potential is presented in Sect. II. There 
we also discuss how to calculate the interaction associated with the convection current. 
We propose in Sect. Ill a simple truncation in the the sum over off-diagonal multipole 
components. We test this approximation in Sect. IV in first-order perturbation theory 
for distant collisions by comparing to the exact results of Ref. Sect. V contains our 
conclusions. 

II. ELECTRO-MAGNETIC INTERACTION WITH FAST CHARGED PARTICLE 

The electro- magnetic interaction of a fast charged particle (with atomic number Zi) and 
a weakly bound proton in a target nucleus has the form 

Km(r, r') = Z,e'{^ _ _I_[p0 + 0p]). (i) 

The charged particle is assumed to move with constant velocity v in the z-direction on the 
straight-line trajectory r' = b + vt, with impact parameter b with respect to the target 
nucleus. The coordinates of the proton in the target nucleus are denoted by r, and p is the 
conjugate momentum operator. The potential is the Lienard-Wiechert potential 

</> = ; ^ (2) 

w |r_L — r^p + 7^(2; — z'y 

where 7 = — (3'^, and (3 = v/c. The second part of the interaction (1), which contains 

the momentum operator p, is due to the convection current. We ignore the magnetization 
current but it could be included by replacing the momentum operator by p + h'V x /x^, 
where /x^ is the magnetic moment of the proton. 

In the two-body breakup of a target nucleus, into a proton and a core nucleus, one 
would actually have to consider the effect of two interactions with the projectile, namely, 
the 'direct' interaction with the proton and the 'recoil' interaction with the core. We will 
not discuss this complication here, since the multipole expansion of the two interactions can 
easily be generated from the formulation given below. 
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A. Multipole expansion 

The multipole expansion of the Lienard-Wiechert potential can be obtained from the 
Fourier representation P 

47r r , e^'»-('-'"') 47r f , e^'^'^'-'"') 1 



(27r)3 J ^ ql + ql/-^^ (27r)3 7 ^ g2 i _ ^2 ^032 

Except for the last factor, this is just the ordinary Coulomb potential. To proceed, we 
introduce the multipole expansion of the last term in Eq. (3) 

^^^=Y.smPM^), (4a) 

1^ ^ ' A even 

where 

^,(/?) = (2A + l)/3-^Q,(/5-^), (46) 

and 

are Legendre polynomials of the second kind; cf. Eq. 8.825 of Ref. or Eq. (8.8.3) of Ref. 
||1CI|| . Explicit expressions and recursion relations for Qx{z) are given in Ref. |TD[. 
We can now insert the plane wave expansion 

Im 

for the two plane waves in Eq. (3) and obtain 

= 47r ^ i'-' Ru>{r,r') A;^,zw(/5) Yil{r) F,w(r'), (5) 

Iml'm' 

where 

2 f°° 1 f°° dq 

Ri,i'{r,r') = - dq ii{qr) ii,{qr') = ^= / — J«+i/2(gr) J;/+i/2(gr'), (5a) 
Ti Jo \Jrr' Jo q 

and 



Am,/w(/3) = {Yi'^m'\y—^^^—^AYi^) = E {Yv^'{q)\Px{0<i)\Yimm. (56) 

i p cos [Uq) ^ g^g^ 



The above expressions were derived for a trajectory with constant velocity in the z- 
direction. The matrix (5b) is therefore diagonal in m. Let us also give a more general result 
which does not refer to any specific coordinate system. This can be done by noticing that 
the angle 9q in Eq. (5b) is actually the angle between q and the velocity v. Thus by inserting 

io(«)*^io( 



Px{0,) = E,D'^,{vyD'^,{q) we obtain 



Aim,i'm'{(3,v)= E 9x{f3) EKo(^^)) {lmXfi\l'm') {1'0X0\10), (5c) 

A ever M 

where the last two Clebsch-Gordan coefficients is the matrix element {Yi>m'\D^Q\Yirn). 

We note that the expansion (5) contains diagonal (i.e. / = I') as well as off-diagonal 
multipole components (i. e. (Im) ^ {I'm'). The off-diagonal components are caused by 
the Lorentz contraction which destroys the spherical symmetry of the potential (f). The 
radial dependence is determined by the integral (5a), and analytic expressions and recursion 
relations are derived in App. A. In numerical applications one can include the finite size of 
the projectile simply by multiplying the integrand in Eq. (5a) by the Fourier transform of 
the charge distribution. However, we find it instructive to show the analytic expressions for 
point-particles. 

The radial dependence of the diagonal components is determined by (see Eq. (A. 5)) 

where r< = min(r, r') and r> = max(r, r') . This is identical to the dependence one has 
in the non-relativistic limit. In fact, the well-known multipole expansion of the Coulomb 
interaction is recovered from Eq. (5), when we insert the non-relativistic limit of Eq. (5b), 

viz. Am,;'m'(/3 = 0) = 5ii>5mm'- 

The off-diagonal multipole contributions to are more complicated. The radial form 
factors have the symmetry property 

Ri,,{r,r')^R,,i{r',r). (7a) 

According to Eq. (5b), we only need expressions for even values of \l — l'\- We show in App. 



A that 

Ri,i+A{r, r') = 0, for r > r' and A = 2, 4, ... (76) 

The exphcit, non-zero expressions can be obtained from Eq. (A. 4). For A = 2 we obtain in 
particular (see Eq. (A. 5)) 

Wr,r') = ^^[l-(^n, for r<r'. (8) 

One can exploit the properties (7a) and (7b) and write expressions for that are valid 
for distant and close collisions, respectively. Thus we obtain for distant collisions (r < r') 

0dist = 47rY;;^(f) Rl,l+^ir,r') Aim,l+Am' Yi+A^rnir'). (9a) 

Imm' A=0,2,. 

For close collisions (r > r') we can use the expression 

0close = E 47rr4(f) Y Rl,l-A{r,r') Aim,l-Am' 11-A,m'(r'). (96) 
Imm' A=0,2,. 

The sum over A is finite for close coUisions, since I — A must be non-negative. This is a very 
nice feature, which makes it feasible to include all terms. For distant collisions, on the other 
hand, one would have to make a truncation in the sum over A. 

B. Convection current interaction 

The contribution to the interaction (1) that originates from the convection current is 

The operator v • V is effectively of dipole nature, as shown in App. B. In practical numerical 
calculations it will operate on an expression of the form r~^f{r)YLM{r)- ^From Eq. (B.4) 
we then obtain 

'fir) 



r fJ, L'M' 

x(i^-Sli'(i'+l)-Mi+l)l). (11) 
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We can also give an explicit expression for the term c • (V0), which appears in second 
version of Eq. (10). Assuming that v points in the z-direction, we obtain from Eq. (3) 

rf0 47r f e^^^''^^') iPqcos{9g) 
^^d^' ~ (2^y ^ g2 l-/32cos2(^,)- 

We can now repeat the derivation that lead to Eq. (5) and obtain 

\^^J IV m 



where 



and 



2 /■°° 

Pii'{r,r') = - dqq ji{qr) ji'{qr'), (12a) 

TT Jo 



Bu'mm = {Yl'm\y^^^%^AYlm) = E 9x{P) (l^'n^l Pa(cOs(^)) J . (126). 
i p cos [t^) ^ 

The last expression is similar to Eq. (5b) but the sum is now over odd values of A. The 
dependence on the radial coordinates can be derived from Eqs. (A. 7-8). We note that matrix 
elements of 4>' can be calculated in a much simpler way in first-order perturbation theory, 
as we shall see in Sect. IV. 



III. SIGNIFICANCE OF RELATIVISTIC CORRECTIONS 

In order to illustrate the significance of relativistic corrections at intermediate energies, 
we show in Fig. 1 the functions gx{P) defined in Eq. (4b). They are seen to decrease rapidly 
as function of A when (f/c)^ < 0.5. This suggests that one would only need a few terms in 
the sum over A in Eq. (9a). 

In the next section we determine the range of velocities where a truncation to A = 0, 2 
is reasonable for distant collision. Thus we will apply the approximation 

~ 21 + 1 

X E[^Wm'(/3) Yirn'ir') - Aw+2,„^'(/5) YwMr') (1 - ir/r'f)]. (13) 

m' 



We assume a straight line trajectory with constant velocity in the z-direction and calculate 
the Coulomb excitation in first-order perturbation theory. We can then test the accuracy of 
the truncation, Eq. (13), by comparing to the exact results from Ref. [|1|. 



IV. TEST IN FIRST-ORDER PERTURBATION THEORY 

In first-order perturbation theory for distant collisions one needs the Fourier integral of 
the multipole fields, 



oo / Am- 

e^-V,^(r,t) = r' F4(f) SiU^), (14a) 



where hu is the excitation energy and Sim{i^) are the so-called orbital integrals. 



For non-relativistic Coulomb excitation and a straight-line trajectory, these integrals are |I| 

Srj^iu) = I ''"I (^)' K^{.hlv). (15) 

^ ^J {I + m)\{l - m)\ ^ 

Below we derive the orbital integrals for the truncated interaction (13). 



A. Lienard-Wiechert potential 

The Fourier transform (14a) of the multipole field (13) can be expressed in terms of 
various combinations of the non-relativistic orbital integrals (15). This is obvious for the 
first term in (13), and the second term is calculated in App. C. A slight complication is the 
very last term in Eq. (13) which contains the factor (r/r')^. It turns out to be a factor of 
{(jjr/vY smaller than the dominant term so let us ignore it. Thus we obtain 



The latter integral is given in Eq. (C.5). 
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B. Convection current interaction 



We also need to calculate the contribution from the convection current, Eq. (10). Here we 
can again make the substitution (B.l). Assuming that the single-particle potential commutes 
with r we can go a step further and replace the momentum operator by Pz = ^ [H, z] , where 
H is the single-particle Hamiltonian. Thus we obtain 

0' = -^([^' ^]<f> + m ^]) = -^([^' ^<f>] + ^Hz - zHcf). (17) 

We note that in first-order perturbation theory one can effectively replace the commutator 
[H, zcf)] by hujz(j). The other term, (pHz — zH(p, has to be considered more carefully. 
Let us write the multipole expansion of (p' as in Eq. (14a), 



Since z in Eq. (17) is a dipole operator, we see that can receive contributions from both 
<Pi-i,m and (pi+i^m- The contribution from the latter is smaller (by a factor of (cur/v)^) so let 
us ignore it. Let us first consider the contribution from the commutator [H, zcj)]. This leads 
to the expression 

SLi[H, #]) - ^ ^.-1,. ^f^. (18) 

For the dipole field we see that the second term in Eq. (17), namely (pooHz — zi?0oo, is 
identical to the first term because 0oo commutes with the Hamiltonian H. Thus we obtain 

S[o ~ -P^ -Soo, S[, = 0. (18a) 
c 

For the quadrupole field we can neglect the second term of Eq. (17) It is identical to zero 
when we insert 0io, whereas inserting leads to a magnetic type transition (oc {l^ + ily))- 
Thus we will use the approximation 



« -f ^ S.,™ (186) 
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C. Comparison to exact results 



In order to illustrate the orbital integrals we obtain from the approximation (13), it is 
useful to plot the reduced values 

Sim = Where iV,^ = - ' ' (19) 

The results we obtain for dipole and quadrupole excitations are shown by the solid symbols 
in Figs. 2 and 3, respectively, as functions of the adiabaticity parameter ^ = ujh/v. Here we 
have chosen the beam velocity v/c = 0.6 

The full relativistic expression for electric transitions which also includes a contribu- 
tion from the convection current, is 

Sf^iuj) = S^Jt/Nim = - GiUf3) (20) 

7 7 

where Gim can be extracted from [^, 

1 

Gio = G20 = G2±2 = —J Gi±i = 1, G2±i = 1 — --. 

7 2 

This result is shown by the solid curves in Figs. 2 and 3. The dashed curves represent the 
non-relativistic result = K^iOy c. f. Eq. (15). It is seen that relativistic effects are 
significant and that the truncated calculation is in very good agreement with the full result. 
If we choose a larger velocity, say v/c = 0.8, some discrepancy starts to occur, so it will be 
necessary to include higher values of A in the sum (9a). 



V. CONCLUSION 

The multipole expansion of the electro-magnetic interaction in relativistic heavy-ion col- 
lisions that we have presented is exact for a straight-line trajectory with constant velocity. 
It has not previously, as far as we know, been used in calculations of the Coulomb excita- 
tion. The formulation is, in particular, applicable to reactions where close collisions play 
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a role, and it is also well suited for calculations of higher-order processes in the Coulomb 
dissociation of halo nuclei at intermediate energies. 

The multipole expansion contains diagonal as well as off-diagonal components, associated 
with the intrinsic coordinates and the coordinates for the relative motion of the colliding 
nuclei. The off-diagonal multipole components are caused by relativistic effects, essentially 
by the Lorentz contraction of the Coulomb field; only the diagonal components survive in 
the non- relativistic limit. 

The number of off-diagonal multipole components is finite for close collisions but it is, 
in principle, infinite for distant collisions. However, it is sufficient to include just a few off- 
diagonal terms at intermediate energies. Thus we find that the diagonal and first off-diagonal 
terms (i. e., A = 0, 2 in Eq. (9a)) give excellent results in first-order perturbation theory for 
distant collisions and for velocities up to f/c ~ 0.6, i. e. up to about 200 MeV/nucleon. 
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where = / + 1/2 and 2m is an even number, according to parity selection of Eq. (5b). 
This integral is given in many textbooks but we shall mainly refer to Abramowitz [|l^. We 
note that it has the symmetry 



VI. APPENDICES 



A. Bessel function integrals 



Here we study the integral in Eq. (5a) 




(Al) 



{a,b) = I^,,^,+2m{b,a) 



{A.2) 



so it is sufficient to determine (A.l) for all values of a and b. 
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Let us first assume that a < b. From Eq. 11.4.34 of Ref. we obtain (witli u = /i + 2m 
and A = 1) 



[a,b) 



1 ,CL, 



r(/i + m) 



2Fi(/x + m, -m; /x + 1; (a/bY). (A.3) 



2'6' r(/i + 1) r(m + 1) 
Tlie liypergeometric function in Eq. (A. 3) is a polynomial of degree m, which can be 
obtained directly for the series Eq. 15.1.1 of Ref. [1^. Inserting into that expression the 
Pochhammer's symbols (see Eq. (6.1.22) of Ref. [1^]) we obtain 



'a,b) 



la 



^-1)'^ T{fi + m + n) 



For m = 0, 1 and 2 we obtain in particular 



(a/6) 



n! 



2n 



(A.4) 



I,A^,b) = j-ia/br, 



1 



{a,b) = -{a/by\l-{a/b) 



W4(a, b) = -{a/br [1 - ia/by\ + 1 - (/. + 3)(a/6) 

Using the symmetry property (A. 2) we see that /^^^ is consistent with Eq. (6). 
When a > b and m is non-zero we find that 

I^i,^l+2m{a,b) = 0, for a>b, m = l,2,... 



(A.5) 



(^.6) 



This result follows from Eq. 11.4.33 of Ref. because that expression will contain the 
factor l/r(— m + 1) which is zero for m = 1, 2, 3, ... 

The functions P^ir{r,r') defined in Eq. (12a) can be obtained from the relati( 



tion 



21' + 1 



{A.7) 



This expression follows directly from the definitions in Eqs. (5a) and (12a) and the 
relation for spherical Bessel functions, zji^i{z) + zji^i{z) = {21 + l)ii{z). For V = I + 1 
find explicitly (c.f. Eq. 6.575 of Ref. §) 



recursion 

we 



1 , r 



Pi,i+i{r,r') = Pi+i^i{r',r) = for r <r' 



(A.8) 



while it is zero when r > r'. The radial dependence can now in general be derived from the 
recursion relation (A.7). 
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B. Current operator 



Here we derive an expression for the operator v ■ V, which appears in the contribution 
from the convection current. It can be expressed in terms of a commutator with the Laplace 
operator A as follows 

vV = ^[A,vr]. (B.l) 

Here we can insert 

V • r = vrcos{e,r) = vrY^{Dlo{v)y Dl,{f), {B.2) 
and express the Laplace operator in spherical coordinates 

where L is the angular momentum operator. Thus we obtain 

V ■ V = vY.{Dl,{v))\Dl,{f){j- + i) - l[LM^io(r)]). (5.4) 



C. Orbital integral 

Here we derive an expression for the integral 



Sn.M = I dt e^- (ai) 



21 + 1 j-oo 

which appears in Eq (16). This integral can be calculated in a simple way by observing that 



hi 4- ^ 1 



where (use Eq. (14) of App. D in Ref. [rT|) 



;(/ + 2 + m)(/ + 2-m) 
1 + 2 

1 ^(/ + l-m)(/ + 2-m) 
' V2 1 + 2 

^Jil + l + m){l + 2 + m) 
'~V2 1 + 2 
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The trajectory we consider falls in the x-z plane, with cos(^) — vt/R and sin(^) = h/R. 
The integral (C.l) can therefore be expressed in terms of the non-relativistic orbital integrals 
(15) as follows 

Sl2jn{^) - 



121 + 5 



21 + 1 



d 



cl v{-i—) s(ilM 



Working out the details this can be written explicitly as 

Al+2+m 



{C.3) 



1 21 + 5 1 



0' 



V V 2Z + 1 1 + 2 ^(Z + 2 + m)!(Z + 2-m)! 

[(2/+3X((/+2-m)X^+i(0+(^+2+m)X^_i(0)-2(/+l)(/+2+m)(i+2-m)X^(0], (C.4) 
where ^ = ujh/v. Or even more compact: 

2 l2l + b 



d+2+m 



V V 2Z + 1 ^(Z + 2 + m)!(Z + 2-m)! 



;-)' 



2? I Q 

{{I + m + 2)+ m')K^{0 - + K^+iiO) 



(C.5) 
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FIGURES 

FIG. 1. The functions gx{v/c), defined in Eq. (4b), axe shown as functions of (v/c)^ for A = 0-6. 

FIG. 2. The reduced values of orbital integrals, defined in Eq. (19), are shown for dipole 
excitations as functions of the adiabaticity parameter ^ = oob/v. The beam velocity is v/c = 0.6. 
The solid curves are the exact results, Eq. (20), and the dashed curves are the non-relativistic 
results. The solid symbols are the results of our approximation discussed in Sect. IV. 



FIG. 3. Similar to Fig. 2 but for quadrupole excitations. 
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